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Dielectric four-port devices play an important role in op- 
tical quantum information processing. Since for causality 
reasons the permittivity is a complex function of frequency, 
dielectrics are typical examples of noisy quantum channels, 
which cannot preserve quantum coherence. To study the ef- 
fects of quantum decoherence, we start from the quantized 
electromagnetic field in an arbitrary Kramers-Kronig dielec- 
tric of given complex permittivity and construct the trans- 
formation relating the output quantum state to the input 
quantum state, without placing restrictions on the frequency. 
We apply the formalism to some typical examples in quan- 
tum communication. In particular we show that for Fock- 
entangled qubits the Bell-basis states |^*) are more robust 
against decoherence than the states I*!?*). 
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I. INTRODUCTION 

Quantum communication schemes widely use dielec- 
tric four-port devices as basic elements for constructing 
optical quantum channels. A typical example of such 
a device is a beam splitter as a basic element not only 
for classical interference experiments but also for imple- 
menting quantum interferences. Another example is an 
optical fiber, which can be regarded as a dielectric four- 
port device that essentially realizes transmission of light 
over longer distances. 

Dielectric matter is commonly described in terms of 
the (spatially varying) permittivity as a complex func- 
tion of frequency, whose real and imaginary parts are 
related to each other by the Kramers-Kronig relations. 
Since the appearance of the imaginary part (responsible 
for absorption and/or amplification) is unavoidably asso- 
ciated with additional noise, dielectric devices are typical 
examples of noisy quantum channels. Using them for gen- 
erating or processing entangled quantum states of light, 
e.g., in quantum teleportation or quantum cryptography, 
the question of quantum decoherence arises. 

In order to study the problem, quantization of the elec- 
tromagnetic field in the presence of dielectric media is 
needed. For absorbing bulk material, a consistent formal- 
ism is given in 0| , using the Hopfield model of a dielectric 
. A method of direct quantization of Maxwell's equa- 
tions with a phenomenologically introduced permittivity 
is given in It replaces the familiar mode decomposi- 
tion of the electromagnetic field with a source-quantity 



representation, expressing the field in terms of the classi- 
cal Green function and the fundamental variables of the 
composed system. The method has the benefit of being 
independent of microscopic models of the medium and 
can be extended to arbitrary inhomogeneous dielectrics 
. All relevant information about the medium are con- 
tained in the permittivity (and the resulting Green func- 
tion), and quantization is performed by the association of 
bosonic quantum excitations with the fundamental vari- 
ables. 

Quantization of the phenomenological Maxwell field is 
especially well suited for deriving the input-output re- 
lations of the field on the basis of the really ob- 
served transmission and absorption coefficients. In par- 
ticular, there is no need to introduce artificial replace- 
ment schemes. Applications to low-order correlations in 
two- photon interference effects have been given |^,|l^,|ll| . 

The formalism has also been extended to amplifying 
media (|Jl^. The resulting input-output relations for 
amplifying beam splitters have been used to compute 
first- and second-order moments of photo counts [ p^ and 
normally ordered Poynting vectors Further, prop- 

agation of squeezed radiation through amplifying or ab- 
sorbing multiport devices has been considered |15| . 

For the study of entanglement, however, knowledge of 
some moments and correlations is not enough. In par- 
ticular, to answer the question as to whether or not a 
bipartite quantum state is separable and to calculate the 
degree of entanglement of a nonseparable state, the com- 
plete information on the state is required in general. 

Recently we have presented closed formulas for calcu- 
lating the output quantum state from the input quan- 
tum state fl^ , using the input -output relations for the 
field at an absorbing four-port device of given complex 
refractive-index profile. In this paper we apply these re- 
sults to study the entanglement properties infiuenced by 
propagation in real dielectrics and extend the theory also 
to amplifying four-port devices. Enlarging the system 
by introducing appropriately chosen auxiliary degrees of 
freedom, we first construct the unitary transformation 
in the enlarged Hilbert space. Taking the trace with re- 
gard to the auxiliary variables, we then obtain the sought 
formulas for the transformation of arbitrary input quan- 
tum states. Finally, we discuss some applications, with 
special emphasis on the dependence of entanglement on 
absorption and amplification. 

The paper is organized as follows. In Sec. |l^ the basic 
equations are reviewed and the general transformation 
formulas are derived. Examples of possible applications 
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are discussed in Sec. ED, and some conclusions are given 
in Sec. fvl. 



II. QUANTUM STATE TRANSFORMATIONS 

A. Basic equations 

Let us briefly review some basic formulas needed for 
the following calculations. For simplicity, we restrict our- 
selves to a quasi-one-dimensional scheme (Fig. The 
action of the dielectric device on the incoming radiation 
is described by means of the characteristic 2x2 trans- 
formation and absorption matrices T(a;) and A(w) re- 
spectively, which are derived in |^ on the basis of the 
quantization scheme in j^. They are given in terms of 
the complex refractive-index profile n(x, lu) of the device. 
Let ai{uj) and bi{uj), i — 1,2, be the amplitude operators 
of the incoming and outgoing damped waves at frequency 
Lo. Taking their spatial arguments at the boundary of the 
device, we may regard them as being effectively bosonic 
operators Further, let gi{uj) be the bosonic operators 
of the device excitations, which play the role of opera- 
tor noise forces associated with absorption or amplifica- 
tion. Introducing the two- vector notation a(aj), b(ci;) and 
g(a;), for the field and device operators respectively, we 
may write the input-output relation for radiation at an 
absorbing or amplifying device in the compact form 



b(cj) = T(cj)a(Lj) + A{uj)d{uj), 



(1) 



where the transformation and absorption matrices satisfy 
the relation 



T(w)T+(w) + aA{uj)A+{uj) = I, 



(2) 



and (7 = +1, d(cj) = g(a;) for absorbing devices and 
cr = — 1, d{uj) — g'iuj) for amplifying devices. The 
above given equations are valid for any chosen frequency. 
Knowing the amplitude operators as functions of fre- 
quency, the full-field operators can be constructed by 
appropriate integration over the frequency in a straight- 
forward way [p| . 



B. Unitary operator transformations 

The operator input-output relation (|l|) contains all the 
information necessary to transform an arbitrary function 
of the input-field operators into the corresponding func- 
tion of the output-field operators. In particular, it en- 
ables one to express arbitrary moments and correlations 
of the outgoing field in terms of those of the incoming field 
and the device excitations, and hence all knowable infor- 
mation about the quantum state of the outgoing field can 
be obtained. Commonly, quantum states are expressed 
in terms of density matrices or phase-space functions - 



representations that are more suited to study a quantum 
state as a whole. 

In order to calculate the density operator of the out- 
going field for both absorbing and amplifying devices, we 
follow the line given in jl^ for absorbing devices. We 
first define the four- vector operators 



a.{uj) = 



a(w) 
d(c.) 



/3H = 



h{u;) 



(3) 



where i{uj) — h{uj) for an absorbing device, and i{uj) — 
W (uj) for an amplifying device, with h(w) being some 
auxiliary bosonic (two-vector) operator. The input- 
output relation (^) can then be extended to the four- 
dimensional transformation 



f3{Lu) = A{Lu)a{u) 



with 



A{uj)JA+{uj) = J, 



J = 



I 

(tI 



(4) 



(5) 



Hence, A(a;)e SU(4) for absorbing devices |T^, and A(w) 
€ SU(2, 2) for amplifying devices (if an overall phase fac- 
tor is included in the input operators). Note, that lossless 
devices, where A{uj) = 0, can be described by SU(2) group 
transformations Since the group SU(4) is com- 

pact, while SU(272) is noncompact, qualitatively differ- 
ent properties of the state transformations are expected 
to occur in these two cases. Introducing the (commuting) 
positive Hermitian matrices 



c(^) = Vtht+H, s{lu) = Vaha+h, (6) 

which, by Eq. (||), obey the relation (uj) + aS^ (uj) =1, 
it is not difficult to generalize the matrix A(ti>) in llq] to 



Ml.) -I ^^^^ ^^^^ 



(7) 



Both the SU(4) and SU(2, 2) group elements can be 
written in exponential form 

A{Lu)^e-'^^^\ ^+{lu)^J^{lu)J, (8) 

and a unitary operator transformation 

= U^a{uj)U (9) 

can be constructed, where 



U — exp 



(10) 



Note that the unitarity of U follows directly from Eq. (H). 

Let the density operator of the input quantum state be 
a functional of a.{(jj) and a.\uj), — gin[Q:((jj), a^(cj)]. 
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The density operator of the quantum state of the outgo- 
ing fields can then be given by 



Tr(-°){gin JA+{uj)Ja{Lj),JA^{Lj)Ja''{uj) }, (11) 



Sout 



where Tr'-'^-* means trace with respect to the device vari- 
ables. It should be pointed out that g^'^l in Eq. ( |ll| ) 
does not depend on the auxiliary variables introduced 
in Eq. (^). The SU(4)-group transformation preserves 
operator ordering and thus for absorbing devices, the s- 
parametrized phase-space functions transform as 



Pout[a(t^);s] = Pin[A+ {uj)a{uj); . 



(12) 



Since the SU(2,2)-group transformation mixes creation 
and annihilation operators, an equation of the type ( |l^ ) is 
not valid for amplifying devices in general. An exception 
is the Wigner function that corresponds to symmetrical 
ordering (s = 0): 



Wout [a{io)] = Win[J A+{uj) Ja{Lu)]. 



(13) 



For amplifying devices, the calculation of the output 
state is rather involved in general. Formulas for Fock- 
state transformation are given in the Appendix. 



III. APPLICATIONS 

As already mentioned, the input-output relation (Q) 
enables one to calculate arbitrary moments and corre- 
lations of the outgoing field. It is worth noting that 
there is no need to introduce fictitious beam splitters 
for modeling the losses. The transmittance and absorp- 
tion matrices in Eq. (^ automatically take account of 
the losses, because they are calculated from Maxwell's 
equations with comple x perm ittivity. To give an exam- 
ple, we compute in Sec. [II A the visibility of interference 
fringes in photon-number detection in a Mach-Zehnder 
interferometer with lossy beam splitters. 

The input-output relation ( pj] ) can advantageously be 
used when knowledge of the transformed quantum state 
as a whole is required. This is typically the case in quan- 
tum communication, which is essentially based on en- 
tangled quantum states. For quantification of entangle- 
ment - a quantum-coherence property that sensitively 
responds to losses - information about the full quantum 
state is needed in general. The entanglement measure we 
use is the quantum relative entropy (the quantum analog 
of the classical KuUback-Leibler entropy) defined by p9| ] 



E{a) — min Tr [a (In — In p)] 

pes 



(14) 



with (7 and S being, respectively, the bipartite quantum 
state under study and the set of all separable quantum 



states. We stress here that the relative entropy is in- 
deed a "good" entanglement measure, because it satisfies 
the necessary conditions that should be required of such 
a measure fl^-^. Note that any proper entanglement 
measure satisfying them would do (the Bures metric be- 
ing another typical example). 

In Sec. Ill B we study the entanglement produced at a 
realistic beam splitter by initially uncorrelated photons, 
and in Sec. [II C we analyze the degradation of entan- 
glement during propagation through lossy media, with 
special emphasis on Bell-type states. Effects associated 



with amplification are addressed in Sec. HID. 



A. Visibility of interference fringes 

To give an example of application of the input-output 
relations (|l|), we consider the visibility of interference 
fringes in a Mach-Zehnder interferometer in Fig. ^ A 
single photon is fed into one input port, the other input 
port being unused. The quantity we are interested in is 
the visibility 



[rik) 



[rik) 



{rik) 



+ (nk) 



(15) 



where (fT.fc)max {{fi-k)min) IS thc maximum (minimum) 
value of the mean photon number in the fcth output chan- 
nel (fc = l,2). 

In order to model the losses in the interferometer arms 
(e.g., nonperfect mirrors or dissipation processes in opti- 
cal fibers connecting the beam splitters BSl and BS2), in 
p2| a fictitious (nonabsorbing) beam splitter is inserted 
into each branch of the interferometer. In practice, how- 
ever, the beam splitters BSl and BS2 are also expected 
to give rise to some losses. Whereas the losses arising 
from the beam splitter BSl may be thought of as being 
included in the replacement scheme considered in [ p2| , 
inclusion in the calculation of the losses arising from the 
beam splitter BSl would require that two additional fic- 
titious beam splitters were inserted between the beam 
splitter BS2 and the detectors. Altogether, a replace- 
ment scheme with four fictitious beam splitters at least 
must be considered in order to model all the losses. 

Application of the input-output relations (^ shows 
that there is no need for such an involved replacement 
scheme. Instead, the proper transmittance and reflection 
coefficients of the beam splitters and mirrors (or fibers) 
can be used to obtain the correct physics, including the 
losses. Applying the input-output relations (|l|) succes- 
sively to the beam splitter BS2, the lossy branches, and 
the beam splitter BSl and assuming the devices are in 
the vacuum state, so that the overall input state is IV'in) 
= |1, 0, 0, 0), it is not difficult to show that 



{h,) = \R,m\w 



|ri|2|T4nT2p 



-^2|i?i||i?2||Ti||T2||T3||T4|cosei, 



(16) 



3 



^2 



|i?inT3|2|r.|2 



21 ^ii?2nT4nrip 



-2|i?i||i?2||Ti||T2||r3||r4| coses. 



(17) 



where 91 = 9 + ipa^ - (fTi + VR2 ~ V'Ta + - <ifT4 and 
92 = 9 + ipR^ - ipTi - 'PR2 + 'PT2 + VT3 - fTi ■ Here and in 
the following, the notation (T;)ii = (T;)22 = Ri = le*"^"' 
and {Ti)i2 = (T,)2i =Ti = \Ti\e''^'^i for the elements of the 
transmittance matrix T/ of the Zth four-port device is 
used [/ = 1, 2, beam splitters BSl and BS2; Z = 3, 4, upper 
(3) and lower (4) branch of the interferometer]. Note 
that for a lossy device 

argi?z-argr, ^ ^2 (18) 

in general. Combining Eqs. ( [T5| ) - (p7|), we easily derive 

|i?i||T3||i?2| , |ri||r4||r2| 



Vi 



V2 



ri||r4||r2| 
(\Ri\\mT2\ 



l^l||T^3||it!2| 
|fi2||Ti||T4| 



V|i?2||Ti||r4| |i?i||T3||r2 



(19) 
(20) 



It is worth noting that Eqs. (|T^) and (|20| ) are valid for the 
really observed reflection and transmission coefflcients Rk 
and Tfc, respectively, with 



\Rk\^ + \Tk\^ <1. 



(21) 



The equality sign would be realized for nonabsorbing de- 
vices. 

Comparing with the formulas for the visibilities de- 
rived in ll^, we observe that they look like Eqs. ([l^) and 
(pO|). However, this resemblance is only formal. In fact, 
all the reflection and transmission coefficients (including 
the phases) introduced in satisfy the relations valid 
for nonabsorbing devices and therefore differ from the 
measured reflection and transmission coefficients that in 
a real experiment determine the fringe visibilities. Even 
if additional fictitious beam splitters were included in the 
model in | |2^ , there would be no unique relation between 
auxiliary and actual parameters in general. 



B. Photon entanglement at a beam splitter 

Superimposing two nonclassically excited modes by a 
lossless beam splitter, one can generate entangled states 
with interesting properties [0 . Two of the simplest ex- 
amples are as follows. Having a single-photon Fock state 
in one input channel of a 50%-50% beam splitter, i.e., |Tp 
= = 1/2, whereas the other input channel is unused, 
the output state is a superposition of states with the pho- 
ton in one of the output channels. If each of the two in- 
coming modes is prepared in a single-photon Fock state, 
then the output state is a superposition of states with 
two photons in one output channel. In either case, the 
output state is a superposition of two states and the max- 
imum entanglement of In 2 (which corresponds to 1 bit) 
is realized. Note that for pure states the entanglement 



measure ( [T^ reduces to the von Neumann entropy of one 
subsystem. When |Tp^|i?p then the output state in the 
latter case is a superposition of three states, because each 
outgoing mode can now contain either zero, one, or two 
photons. The maximum entanglement of a three-state 
system is In 3, which is realized if \T\^ = 1/2 x (1 ± l/\/3) 
(see Fig. ^). Hence, a non-50%-50% beam splitter can 
produce stronger entanglement than a 50%-50% beam 
splitter which suppresses one possible outcome owing to 
interference. With regard to entanglement, this interfer- 
ence effect is thus destructive. 

Let us now raise the question of the amount of en- 
tanglement achievable in case of a realistic beam split- 
ter - a question that may be important for the quality 
of quantum communication by means of entangled pho- 
tonic states obtained by available devices. The question 
can be answered by applying the input-output relation 
( prj ) and calculating the output state of the interfering 
modes obtained by an absorbing beam splitter. To give 
an example, let us study the entanglement produced by 
a dielectric plate of permittivity 



e{uj) = 1 + 



1 



1 — (uj/ujo)^ — 2ijLu/ujQ 



(22) 



{eg = 1.5) and thickness d — 2c/wo for the case where ei- 
ther one or each of the two incoming modes is prepared 
in a single-photon Fock state. The squares of the abso- 
lute values of the calculated reflection, transmission, and 
absorption coefficients as functions of frequency |^ are 
shown in Fig. ^ for 7 = 0.001. When the device is not 
excited, then the overall input state is either |1, 0,0,0) 
or 1 1, 1,0,0) for the two cases under consideration. The 
resulting mixed states of the outgoing modes are calcu- 
lated in Here, we have calculated the amount of 
entanglement of the states using the definition (p^. 

Results are plotted in Figs, [s] and || for 7 = 0.001, and in 
Fig. ^ for 7 = 0.01. For comparison, the figures also show 
the mutual information Ic = Si + S2 — S'12, where Si and 
S2 are the von Neumann entropies of the outgoing modes 
1 and 2, respectively, and S'12 is the entropy of the com- 
posite two-mode system. Cbviously, the mutual informa- 
tion may be regarded as a measure of the total amount 
of correlation contained in the states. In regions where 
the absorption is sufficiently weak, the output state is 
almost pure, and thus S'12 ~ and E{a) ~ Si. Hence, the 
two curves in Figs. || and ^ differ there only by a factor 
approximately equal to two. With increasing absorption 
the two curves cannot be related to each other by simple 
scaling, as it can be seen from Fig. 0. In particular, the 
maximally achievable amount of entanglement of about 
0.4 is much less than In 2 achievable with a lossless device. 

From Figs. |^ and ^ strong reduction of entanglement 
is observed in the resonance region. Here reflection and 
absorption are strongest, so that the two modes are only 
weakly mixed and absorption prevents the device from 
creating quantum coherence. As expected, substantial 
entanglement is observed in regions where the absorption 
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is weak and jTp and nearly satisfy the condition of 
maximum entanglement. In Fig. |^ this is the case for 
LojujQ Ki 1.25 where the value of entanglement becomes 
close to the maximally achievable value of In 2. In Fig. ^ 
the value of entanglement becomes close to the maximally 
achievable value of In 3 at u/ujo « 1.18 and uj/luq « 1.33. 
The relative minimum in Fig. |^ at uj/ujo~ 1.25 indicates 
the effect of destructive interference mentioned above. 

The results show that entanglement sensitively de- 
pends on the optical properties of the material used for 
manufacturing the optical device. They demonstrate the 
importance of optimizing the frequency regime of quan- 
tum communication schemes with given devices. 



C. Entangled-state transmission through a lossy 
channel 

1. Bell-type basis states l^*^) 

Let US now turn to the question of entanglement degra- 
dation during the propagation through dielectric matter 
such as an optical fiber. For this purpose, we consider 
two modes each of which propagates through a dielectric 
medium of complex permittivity. Assuming the incoming 
modes are prepared in a maximally entangled Bell-type 
state 



1 



^(|On)±|nO)), 



(23) 



we apply Eq. (|ll|) and calculate the quantum state of the 
two outgoing modes. After some algebra we derive 



fout 



ri|"=(i-|Ti| 



E 

k=0 
k=0 

+k{\T,\'- + \T,\'") |vI':±)(vI/'„ 



|A:0)(A:0| 



n 
k 

|2n 



where 



I* 



IT. 



-1/2 



^—k 



|0fc)(0fc| 



-,'±l 



(24) 



T^\nQ)±T^\Qn)). (25) 



Note that when setting n — the transformation of the 
ordinary Bell basis states l^-^^) = are obtained. In 
what follows we assume that the transmission coefficients 
Tfe (fc = 1, 2) are given by 



(26) 



with nfe(w) — y/ek{oj) = r]k{uj) + ink{'-^) and Ik being the 
complex refractive indexes of the media and the propa- 
gation lengths, respectively. According to the Lambert- 
Beer law, \Tk\ decreases exponentially with the length of 
propagation: IT^I = exp(-Zfc/Lfc), Lk = c/{u>Kk)- In spe- 
cial cases when one mode propagates through vacuum. 



n{u}) = 1, the corresponding transmission coefRcient, by 
Eq. (^), is just a phase factor. 

For a first insight into the behavior of the transmitted 
quantum state it may be instructive to look at the overlap 
of the output state with the input state, which is 

(*,tl^iu^l*,t> = 1 (|Tip"+|r2|2"+T*"r2"+ri"Tr) . 

(27) 

We see that the characteristic length of degradation of 
the overlap (fidelity) is not given by but by the shorter 
length Lk/{2n). Hence, the overlap rapidly approaches 
zero with increasing number of photons even for weak 
damping of the intensity or related (classical) quantities. 

As already mentioned, a proper measure of entangle- 
ment is the quantum relative entropy defined by Eq. (M). 
In order to estimate an upper bound, we employ the con- 
vexity property | |23[ | 

E[\ai + (1 - A)o-2] < \E{ai) + (1 - \)E{a2) ■ (28) 
From Eq. ( p^ ) it is seen that gl^J has the form 

gill ^ \a, + {1 - X)a2, (29) 

where cti is a separable state [E{(Ti) =0] and 

<^2 = K^)«±| (30) 

is a pure state, the entanglement of which is simply given 
by the entropy of one of the two modes. We thus find 



B^-[{\T,f'' + \T2\^'')\n{\T,\'^ + \T2 



2n\ 



- |Ti|2"ln|ri|2" - |r2p"ln|T2p"] 
In particular when Ti = T2 = T, then 

EigZi) < |r|2"ln2 = e"2"'/^ln2. 



(31) 



(32) 



(33) 



i.e., the characteristic length of entanglement degrada- 
tion decreases as l/{2ri) at least. The result ( ^3|) reveals 
that with an increasing number of photons the quan- 
tum interference relevant for entanglement exponentially 
decreases at least. Such a behavior is typical of quan- 
tum decoherence phenomena and is not restricted to Fock 
states. 

It should be mentioned that for a pair of spin-^ par- 
ties a decomposition of the density matrix into a sep- 
arable part and a single pure state is always possible 
P^ . Moreover, there exists a unique maximal A such 
that the inequality (|2^) reduces to an equality and thus 
(1 — X)E{a2) becomes a measure of entanglement. How- 
ever, for larger dimensions of the Hilbert space we are 
left with the general inequality 
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Examples of entanglement degradation [calculated on 
the basis of Eq. for singlet states with one photon, 
), and two photons, are shown in Fig. ^ for the 

case where the two modes propagate in equal media over 
equal distances. We observe that for the state l^*^) the 
upper bound e^'"/^ In 2 defined by the inequality (|3^ ) is 
a very good approximation to the entanglement at propa- 
gation length I. In contrast, for the state l^'f ) the actual 
values of entanglement are typically smaller than it might 
be expected from the upper bound e~^'/^ln2. Since for 
n > 2 the upper bound e~^"'/^ In 2 is always smaller than 
the entanglement observed for the state l^*^) (at least for 
< 1<L), we leave with the result that the two-photon 
singlet state |\E'^) is the most robust one within the class 
of states \^^). 



2. Bell-type basis states 
The Bell-type states 



1 



V2 



(|00) ± \nn)) 



(34) 



can be obtained from the somewhat more general class 
of states 



1 



■.{\00) + q\nn)) 



(35) 



for <7 = ±1. Obviously, for n—1 and small values of q the 
state 1$^) approximates a two- mode squeezed vacuum 

exp[C(aia2 - 44)] 1 00) = V'T^^F ^ | m"^> (36) 



Again, from the convexity argument, Eq. (|2§|), an upper 
bound of the entanglement can be derived 



B 



1 



1 



191- 



[(1 



ln(l + |g'| 



/\2\ 



j'plnlgf] 



(39) 



{q' = qTiT2) ■ In particular, for small values of q' we find 
by expansion that 



V|2 



1 + 



ri-ln|g' 



l\2\ 



oiwn, (40) 



.112 



which shows that the entanglement decreases as \q' 

|9|2|Ti|2»|T2|2". 

It is also instructive to compare the entanglement 
degradation of the states 1$^) with that of the states 
I'^n)- Similar to the states l^"^), within the class of 
states |<i>^) the state |$^) is most robust against en- 
tanglement degradation. Obviously, the probability of 
finding n photons in one channel decreases as |Ti|" 
for the states l^*^) but decreases as |rir2|" for the 
states 1$^). The entanglement degradation of the states 
I'l'^) is therefore expected to be less than that of the 
states From Eqs. (||) and (||) it follows that 

(|ri| = |r2| = |r|«i) 



S(|$±)) _ |T|2" (1 -ln|r|4") 
B(|*±)) ~ 21n2 



(41) 



The numerical results (see Fig. ^ indeed show that the 
states 1 4'^) are more robust against entanglement degra- 
dation that the states 



3. Medium with BIT characteristics 



{q — tanhC, C, real) used in quantum teleportation with 
continuous variables . It is not difficult to prove that 
the entanglement of is 



Ei\K){'S>l\)^lnil + \qn 



1 



(37) 



which for I?] = 1 attains the maximum value of In 2. 

Let again consider two modes propagating through di- 
electric matter and assume that the incoming modes are 
now prepared in a state |$^). We again apply Eq. ( |ll| ) 
and calculate the quantum state of two modes. The re- 
sult reads 



^(F) 
fout 



1 + 

X (1 - iTi 



E 

fcl,fc2=0 

2\ n — ki 



n \ n 
ki) W 



|y^|2/ci|y^|2fe. 



1" 



l2\ n— fc2 



|/Ci/C2)(fcifc2 



-|ri|2"|T2p"|nn)(nn| 



X [|00) + qTl^T^\nn)] [(00| + {qT^T^T {nn\] . (38) 



Media having a electromagnetically induced trans- 
parency dispersion characteristics have been of increasing 
interest [^,0 . They may offer the possibility of realiz- 
ing optical quantum gates, because the group velocity re- 
duction is extremely large such that there will be plenty 
of time to manipulate a quantum state intermediately 
stored in the medium |^^. The susceptibility of such a 
medium can be given by 



x{5) = 



^7i(*7o - 5) 



+ 7^70 - 5{A -5)+ i[(5(7_L + 70) + A70] ' 

(42) 



with 17 being the Rabi frequency if the driving field, 7_l 
the transverse relaxation rate of the probe transition, A 
the one-photon detuning, 71 the radiation relaxation rate 
of the probe transition, 70 the decay rate of the ground- 
state coherence, and 5 the two-photon detuning (for de- 
tails, see |2^,||l). 

We have calculated the degradation of entanglement 
for the case where two modes that are initially pre- 
pared in a Bell- type state |^^), Eq. (E3h, can propagate 
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through media of that type. Figure |o| shows results ob- 
tained for ordinary Bell states l^**). In the figure, the 
two-photon detuning is varied in a small frequency re- 
gion around some optical frequency ujq. The two-peak 
structure of the absorption coefRcient [imaginary part of 
the square root of the susceptibility ([42|)] essentially de- 
termines the amount of entanglement that can be trans- 
mitted. It is seen that the initial entanglement of In 2 is 
(approximately) preserved for zero two-photon detuning, 
and the degradation of entanglement is almost abrupt for 
nonzero two-photon detuning. Hence, control of entan- 
glement requires fine tuning. 



D. Entanglement transformation at amplifying 
devices 

From Sec. ^ wc know that quantum-state transforma- 
tion at amplifying four-port devices is connected with 
SU(2,2) group transformations. For each frequency, the 
transformation corresponds to the action of a four-mode 
squeezing operator, where the destruction (creation) op- 
erators of the field modes are mixed with the creation 
(destruction) operators of the device excitations. Trac- 
ing with regard to the device variables then yields the 
(two-mode) output state of the field, as is shown in the 
Appendix for the case where the (two-mode) input state 
of the field is a Fock state and the device is in the ground 
state. 

If the input field is prepared in an entangled state, 
amplification is expected to destroy the entanglement. 
Although all necessary formulas are available, the cal- 
culation of the quantum relative entropy is an effort. 
The number of the (real) parameters specifying an ar- 
bitrary separable density matrix increases dramatically 
with the dimension of the Hilbert space of the subsys- 
tems involved. In fact, it is easy to see that this number 
is [47V* {N-l) + N'^-l], with N being the Hilbert-space 
dimension of the subsystems (here, both subsystems are 
assumed to have equal dimensions). Hence, when there is 
notable amplification, then the number of Fock states to 
be taken into account for sufficient numerical accuracy 
drastically increases. In contrast to absorbing media, 
where the dimension of the Hilbert space of the relevant 
modes is bounded by the total number of input photons, 
such a bound does not exist for amplifying media. 

Nevertheless, for entangled Gaussian states an upper 
bound of the gain can be determined such that the am- 
plified system is still not separable. Let us consider, e.g., 
the two-mode squeezed vacuum (^6|) and assume that the 
two modes travel through amplifying devices at zero tem- 
perature. The Wigner function of the two-mode squeezed 
vacuum is a Gaussian 



W{$) = (47rVdet f) ^ exp|-i^'^F-i$} 



(43) 



Here, ^ is a four- vector whose elements are qi,pi,q2,p2, 
and V is the 4x4 variance matrix 



V = 



X z 
Y 



(44) 



The variance matrix ( [44j) can be written in the form 



V = 



/ c/2 -s/2 

c/2 s/2 

-s/2 c/2 

V s/2 c/2 



(45) 



[c=cosh2C, s = sinh2C]. Using the input-output relations 
(|l|) for amplifying devices, we can easily transform the 
input-state variance matrix ( ^ ) to obtain the output- 
state variance matrix 



V 



fx Zii Zi2 

X Z21 Z22 

Zu Z21 y 

V Z12 Z22 y 



(46) 



where 



lc\T,\' + l\R,\^ + l{\T,\' + \R,\^-l), (47) 



Zii = 

Z\2 



-Z22 = 

- Z2\ = 



\T2? + \R2? - 1) , (48) 



-isRe(rir2), 

isIm(TiT2). 



(49) 
(50) 



Let us consider equal devices, so that Ti = T2 = T and 
i?i = i?2 = R. The Peres-Horodecki criterion [E9| 



det X det Y 



I det Z|)^ - Tr (XJZJYJZ'^j) 



> i(detX + dctY) 



1 
-1 



then tells us that for 



2 



-2|CI 



(51) 



(52) 



the boundary between separability and nonseparability is 
reached. In particular for zero reflection (i? = 0), Eq. ( [5^ ) 
reveals that the upper limit of the gain 17 = |r|^ — 1 > for 
which nonseparability changes to separability is simply 
given by the squeezing parameter \q\, 



g^\q\= tanhlCl ■ 



(53) 



An obvious consequence of Eq. (^3|) is that entanglement 
cannot be produced from the vacuum by amplification. 
Since for the vacuum the squeezing parameter has to be 
set equal to zero, (? = 0, from Eq. (m) it follows that any 
nonvanishing gain g must necessarily lead to a separable 
state. 
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IV. SUMMARY AND CONCLUSIONS 

We have studied the problem of quantum-state trans- 
formation at absorbing and amplifying dielectric four- 
port devices, without making use of any replacement 
schemes. We instead express the input-output relations 
in terms of the actually observed quantities as obtained 
from the quantized Maxwell field in the presence of ar- 
bitrary causal (linear) media. After extending the basic 
formulas recently developed for absorbing media to am- 
plifying media, we have applied the theory to some prob- 
lems typically considered in quantum information pro- 
cessing. 

In particular, we have considered both the amount 
of entanglement that is realized when nonclassical light 
is combined through a lossy beam splitter and the en- 
tanglement degradation when entangled light propagates 
through lossy media. We have based our analysis on the 
quantum relative entropy as a measure of entanglement. 
The calculation of the entanglement of a mixed quantum 
state typically observed for absorbing media needs com- 
paring the state with all separable states in order to find 
that separable state which is closest to the state under 
consideration. Since the effort drastically increases with 
the dimension of the Hilbert space, we have restricted 
our attention to low-dimensional quantum states in the 
numerical calculation. 

The numerical results show that the Bell-type states 
1 5'^), Eq. (p3|), are more robust against decoherence than 
the states [^), Eq. (|^ (n = 1,2). The estimation of 
an upper bound of entanglement for arbitrary number 
n of photons in each of the two entangled modes shows 
that with increasing n the characteristic length of entan- 
glement degradation decreases as L/n at least, where L 
is the absorption length according to the Lambert-Beer 
law. 

So far we have considered either purely absorbing or 
purely amplifying media. In practice, the two effects can 
occur simultaneously. Essentially, there are two ways to 
deal with this problem. One way is to treat amplifiers 
with absorption as cascading amplifying and absorbing 
devices. Another way is to go back to the underlying 
quantized Maxwell equations with the aim to develop a 
more specific approach to the problem. 
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APPENDIX A: DENSITY MATRIX FOR 
GAUSSIAN WIGNER FUNCTION 



As mentioned in Sec. II B in the case of amplifying me- 



dia only symmetric operator ordering is preserved, and 
hence the Wigner function (|l^) is suited to the state de- 
scription. For the sake of transparency we will restrict 
our attention to a single-frequency component [i.e., a 
(quasi-)monochromatic field in a sufficiently small fre- 
quency interval Auj ||l6|] ). The extension to a multifre- 
quency field is straightforward. When the input field 
is prepared in a Fock state \p, q) and the device in the 
ground state |0,0), so that the overall input state is 
\p, q, 0, 0), then the input Wigner function reads as 



xip(4|ai|2)i,(4|a2p)e-2(l'^il^+l«^n 
with Ln{x) being the Laguerre polynomial 



n \ X 
n — ml ml 



(Al) 



(A2) 



We now apply Eq. (g_3|) , making the substitutions accord- 
ing 



a-> T+a- T+C^Sg*, 

g^-A^a*+A^ [S^]"'c^g, 
g* ^ -A+a + A+S^^Cg*. 



(A3) 
(A4) 

(A5) 
(A6) 



Finally, we integrate over the device variables gi to obtain 
the Wigner function of the outgoing field. Introducing 
the matrix Kn' — Sn' ki and employing the formula 



dk 



fc=0 



(A7) 



we derive 



h\ 



l\ 



gh Ql 



exp{-2(a*)^[N - B^(D^)-iB*]a} 



detD 



where the abbreviations 



, (A8) 



fci=fe2=0 



N = 2TT+-I-2TKT+, (A9) 
B = 2S^C^ - 2S*C* -^TRT"^, (AlO) 
D = 2T*T'^ - I - 2S*C* "1t*KT^C^"1S^. (All) 



have been used. 
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In order to calculate from the Wigner function the den- 
sity operator, we make use of the relation 



gZl = TT^ / d'^W^^t' (a, a*) 5(a - a), (A12) 



where 



^(a- a) = — / d^bi:>(b)e'^ 



(A13) 



with D{h) being the two-mode coherent displacement op- 
erator. For notational convenience we introduce the ab- 
breviation notation 



h=0 1=0 



h\ 



Qh Ql 



V+1 



{■■■} 



/! \l 



(A14) 



Substitution of Eq. (|A8|) into Eq. ( |A12D yields 
4 



TT^ det D 



where M = N - B^(D^)-^B+. 

Using the Fock-state representation of the (single- 
mode) coherent displacement operator ]3C[|, 



{m\D{b)\n) ^ 6'"-«e-l''l'/2L('"-")(|fe|2) (A16) 



[L^{x), associated Laguerre polynomial], we can calcu- 
late the density matrix in the Fock basis. Performing the 
a integrals in Eq. ( |A15| ), we derive 



( F) 

(mi,m2|^out |7ii,n2) = V 



1 



7r2 det DM 



E 



ni ,n2 ,mi .ni2 



ni\n2l 
mi!TO2! 



ridri r2dr2 di^i d(p2 



rr 1 



My: 



detM 



1 n / A/ii \ |Afl2| 

X^(m.-«.)(^2)\ (A17) 



where we have used the notation bi — r^e*'''*, and M12 = 
|Mi2|e*®. Recalling the definition of the modified Bessel 
functions, we perform the angular integrals to obtain 



{mi,m2\gi^l\ni,n2) 



V 



1 



det DM 



E 



rj-i ,712 ,'>^'>-l:'>"l"l'2 



ni!n2! 
mi!m2! 



mi— ni+m2— 712,0 



X J dxi dx2 exp 




1 + 



A/2' 



1 



dctM 



detM 

M12 



detM 



X X 



(mi-ni)/2^(m2-n2)/2^(„ii-ni) 



x47-"^)(:e2) 



(A18) 



(xi = T'f ). The X2 integral is performed by means of the 
formula (2.19.12.6) in plf|, which gives (for TO2 > ^2) 



(mi,m2|eiutVi>'^2) 



V 



— y 

D ^ 



det D V mi\m2\ 

ni ,712 ,mi ,m2 



ni!n2! 



(M*2 



m 1 — n 1 + 2 — 2 , 



exp(-2a+Ma + a^b* -b^a*)l 1, (A15) 7 



(Mil - detM)"2 
(Mil + detM)'"2+i 



x / dxi exp 



1^ A I-I-M22 
'2^U Mil + detM 



x4™i-"i)(xi)47-"^) 



Mi\ - (detM) 2 



Xl 



(A19) 



Finally, the xi integral is performed by expanding the 
associated Laguerre polynomials into power series 
The result is 

(TOi,TO2|^i^tVl:"2) = 



V 



— y 

D ^ 



det D ^ V wi!m2! 

ni ,712 , mi ,7772 



ni!n2! 



X (^mi — 



mi— 7ii+m2— 712,0 



(Mil -detM)"2 
(Mil -HdetM)™2+i 



X (M*2)"'""' 



)712 t 
X ^ C / 7712 
Qfc+1 
fe=0 



112- k 



X 2-Fi ( fc + 1, — ni; TOi — rii + 1, - 

a 



where 



1 + Mil + M22 + det M 
2(Mii -hdetM) 
IM12P 



(detM)2-M2^ 



(A20) 

(A21) 
(A22) 



Integrating Eq. (A_8) over the phase space of one mode 
of the outgoing field yields the Wigner function of the 
quantum state of the other mode 
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h=0 1=0 



hi 



d'^ a' detE 
^ dkl^ Eu det D ^ 



■2|ai|VB. 



(A23) 



fcl=fc2=0 



(E = M This Wigner function is equivalent to the 
density matrix in the Fock basis 



= E 

n=0 



EE 

h=0 l=Q 



_-\\h+l+p+q 



(-1) 



hW. 

Eii - 



d^+^ detE 2 



|n)(n|. (A24) 

ki=k2=0 
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FIG. 1. Quasi one-dimensional geometry of the device 
with definition of field input and output operators. 




FIG. 2. Mach-Zehnder interferometer with dispersive and 
absorbing beam splitters BSl and BS2 and a phase shift G. 

The mirrors Ml and M2 as well as the branches between BSl 
and BS2 are assumed to be lossy with transmittance T3 (upper 
branch) and T4 (lower branch), respectively. 
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FIG. 3. Entanglement created at a lossless beam split- 
ter with both input modes excited in an one-photon Fock 
state. The entanglement is shown in exponential scaling vs. 
transmittance. Maximal entanglement In 3 is created at a 
20%/80% beam splitter whereas for a symmetric device it 
drops to In 2 due to destructive interference. 
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FIG. 6. Frequency dependence of the entanglement mea- 
sure E{a) (dashed line) and the total correlation Ic (solid 
line) for a state |1, 1, 0, 0) impinging on a beam splitter with 
7/tJo = 0.001 in Eq. (03) and the other parameters given in 
the text. 
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FIG. 4. The reflection coefficient |_Rp (fuU line), the trans- 
mission coefficient [Tp (dashed line), and the absorption coef- 
ficient (1 — |_R|^ — \Tf) (dotted line) of a dielectric plate are 
shown as functions of frequency ui for es = 1.5 and 'y/ujo = 
0.001 in Eq. (|22|), and the plate thickness 2c/ujo- 




0.5 1 1.5 

FIG. 5. Frequency dependence of the entanglement mea- 
sure E{a) (dashed line) and the total correlation Ic (solid 
line) for a state |1, 0, 0, 0) impinging on a beam splitter with 
7/tJo =0.001 in Eq. (G2) and the other parameters given in 
the text. 
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FIG. 7. Frequency dependence of the entanglement mea- 
sure E{a) (dashed line) and the total correlation 7c (solid line) 
for a state |1, 0, 0, 0) impinging on a beam splitter with 'y/cuo 
= 0.01 in Eq. ( p2|) and the other parameters given in the text. 
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FIG. 8. Entanglement degradation of a singlet state 
l*n)(*nl [Eq. (|2|)] with one photon (n = 1 full curve) and 
two photons (n = 2 dashed curve) after transmission through 
absorbing channels of equal transmittance. 
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FIG. 9. Comparison of entanglement degradation of 
one-photon Bell basis states I"!"*) (full curve) and 
(dashed curve). 
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FIG. 10. Entanglement degradation of a singlet state 
[Eq. ( p3[ ) with n — 1] after transmission of one 
subsystem (dashed curve) or both subsystems (full curve) 
through a medium with susceptibility (f42|). 
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